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For a free scalar boson field and for U(l) gauge theory finite volume (infrared) and other correc- 
tions to the energy-momentum dispersion in the lattice regularization are investigated calculating 
energy eigenstates from the fall off behavior of two-point correlation functions. For small lattices 
the squared dispersion energy defined by i?J is = Ei — E 2 — 4^\_ J 1 sin(fci/2) 2 is in both cases neg- 
ative (d is the Euclidean space-time dimension and E% the energy of momentum k eigenstates). 
With Eo = observation of E 2 iB = has been an accepted method to demonstrate the existence 
of a massless photon in a 4D lattice gauge theory, which we supplement by a study of its infrared 
corrections. A surprise from the free field is that the corrections become larger when the corre- 
lation length £ = m _1 decreases as follows from the derived equation E% — cosh _1 (l + 
M 2 S = M 2 +±YnZl sm{ki/2) 2 , where M is the mass parameter of the free scalar lattice action. 

PACS numbers: PACS: 11.15.Ha 
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I. INTRODUCTION 

The propagator of free fields in the lattice regular- 
ization suggests that the continuum Euclidean energy- 
momentum relation 



becomes replaced by 



eI 



d-l 



El = 



^ + 4^sin(fc i /2) 5 



(1) 



(2) 



Here d is the space-time dimension and E? the energy 

found from correlation functions of momentum k eigen- 
states. Lattice and continuum quantities (the latter here 
with overlines) are related by 



E^ = — and ki = ^ 
k a a 



where a is the lattice spacing. See, for instance, the 
textbook by Rothe [1]. A divergent correlation length 
£ = M^ 1 allows for a quantum continuum limit a — > 0. 

To exhibit corrections to @ we define a dispersion 
energy E dis by 



d-l 



El(N) 



E 2 (N)-Aj2Mh(N)/2}\ (4) 



where N refers to the spatial size of a periodic N d ~ x N d 
lattice. In the following our equations for the free scalar 
field hold for general d = 1, 2, . . . , while for U(l) we per- 
form Markov chain Monte Carlo (MCMC) calculations 
only in 4D. We estimate energy eigenvalues Et through 
the usual cosh-type fits to correlations of operators which 
are in momentum k eigenstates 



Cdn d ) 



,-E s (.N d -n d ) 



(5) 



with parameters and E?. On our lattices rid takes the 
values 0, 1, ... , Nd — 1 and the momenta are discretized 
by 



2tt 

h - o,— , 



o, 



2tt 



2n(N - 1) 
N ' 

2TT(N d - 1) 



1, (6) 
(7) 



We illustrate k ^ with the lowest non-zero momentum 



(fc,0,0) with k = 2ir/N 



(8) 



for which we denote the energy eigenvalue by E^ = E^ 
and the correlation function by Ck{n d )- 

In the next section we derive exact results for the free 
scalar field. In section ITO1 we turn to U(l) lattice gauge 
theory (LGT) in Wilson's regularization Q and rely on 
MCMC investigations. We set E%(N) = and use (E^) 2 



(3) in Eq. fl2J|. Already in the early days of LGT 







(9) 



was estimated in the Coulomb phase for the photon dis- 
persion energy on small 4 3 8 lattices within the numer- 
ical precision available at that time 0. As the direct 
calculation of the zero-momentum photon eigenstate Eq 
fails, the estimate via © became a practical method for 
identifying Coulomb phases in Higgs-type models on the 
lattice. For examples see 0, [f| . 

In Ref. @ the result (|9]) was again consistent with the 
data for simulations in the U(l) Coulomb phase. This 
time from MCMC calculations on much larger 16 4 lat- 
tices. However, in Ref. Q 



disf 



-0.2102 (17) 



(10) 



was reported for a 4 3 16 lattice at (3 = 1.1, where the 
error bar is given in parenthesis and applies to the last 
digits of the estimate (for the Wilson action j3 > j3 c with 



2 



f3 c = 1.0111331 (21) [|| is in the Coulomb phase). In the 
previous work [3j, |6[ the finite size effect was apparently 
each time swallowed by the statistical error (in be- 
cause small systems were not simulated). In section Hill 
we present a systematic study of this finite size behavior 
for which some of the material is taken from the under- 
graduate honors thesis [9(. 

A brief summary is given in the final section HVl 



II. FREE SCALAR FIELD 

Following the standard approach, e.g. [lj, the lattice 
rcgularization for the action of the 4D free scalar boson 
field reads (n and m will be integer four vectors) 



Sr 



^2 ^ <t>n 4>n+p, 
n f_L—l 



M 2 



., (ii) 



where /t is the unit vector in fi direction. Next, we con- 
tinue with a general dimension d and derive analytical 
expressions for the two-point correlation functions on fi- 
nite 7V d_1 Nd lattices with periodic boundary conditions. 
Subsequently, we compare them to the cosh mass fit func- 
tion ©. 

For an infinite lattice the two-point correlation func- 
tion of the scalar field cj> n described by the lattice action 
(fTTjl is derived in [l[ to be 



,J n / oo 



+7r d d k' 



exp[i k' ■ (n — m)] 



w (2^ m 2 +4E d s . n 2 (fc , /2) 



,(12) 



where the integration is over each of the d vector com- 
ponents (d = 4 in lj). Here we are interested in the 
lattice size corrections to this equation. For this purpose 
we replace the integral representation of the Kronecker 
delta 

/+ 7r (jdk 1 
——7 cxp[i k' ■ (n- m)\ (13) 
-7T (27T) d 

used in [l[ by the one for a N d ^ 1 Nd lattice 



N d-x Nd J2 ex P\ ik ' '-(n-m)}. (14) 



Here the summation vectors are given by Eq. ^ and ([7]) . 
To match (|T3|) in the limit N — > oo, Nd — > oo each fc M 
value can be shifted by — ir. Following the logic of [l[ the 
two-point correlation function becomes 



1 



exp[i k' ■ (n — m) 



^- 1 ^^M2 + 4E^ 1 sin 2 (fc;/2) 



(15) 



In this paper we are interested in correlation be- 
tween operators of a definite spacelike momentum k — 



(fci, . . . , fcd-i) for which one expects the functional form 
(O. We define these operators by 



(16) 



where is a normalization constant. We want to calcu- 
late the correlation function 



^^exp ~ik-(n-m) {<p n (j> m ) 



Using translation invariance on a periodic lattice this 
equation simplifies considerably. First, we note that 



^ exp -ik- (n - rn) {(j) n (j> m ) 



E 



exp 



-ik ■ n 



6,mJ 



holds for all m. Besides we have 

( ( t ) n,nd ( l ) 0,m d ) ~ ( ( ( ) fl,n d —md ( f ) 3fi) 

so that it is sufficient to consider 



E ex p 



-ik ■ n 



for which we carry out the n and then k' summations to 
obtain 



H exp(ik' d n d ) 



(17) 



N d j£ AP + 4Ef=i sin 2 (fc l /2) + 4sin 2 (fc^/2) 

The normalization constant 6g can be chosen to ensure 
C^(0) = 1, which fixes also in (TITo)) . From ([T7]) it is ob- 
vious that the momentum k correlation function Cg(n^) 
simply agrees with the momentum zero correlation func- 
tion Co (rid) at the effective mass 



M. 



cfr 



M% +4 2jsin 2 (fc,/2) 



(18) 



It is for the free field sufficient to investigate Co(rid) as 
function of its mass parameter M. Using neighboring 
distances for the correlation function arguments and the 
functional form ([S]), values for Eq are determined by the 
equation 



Co(n d ) 
C (n d + 1 



(19) 



exp[-E n d ] + exp[-g (N d - n d )} 
expf-^o (n d + 1)] + exp[-E {N d -n d - 1)] ' 

Inserting (fT?) for C (rid) and C (n d + 1), the values for 
Eq turn by numerical inspection out to be independent 
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FIG. 1: Lattice size dependence of the dispersion energy (|4|) 
for the free scalar field. 



of rid and Nd for Nd = 2, 3, 4, ... and n<j = 0, . . . Nd — 1 
as one may expect for a massive free scalar field. In 
particular this includes Nd — 2 and nd — for which we 
find 



C o (0) M- 2 + (M 2 + 4)- 1 1 + exp(-2 Eq) 



C (l) M- 2 - (M 2 + 4)- 1 2 exp(-So) 
With a little algebra this simplifies to 

§^ = l + Af 2 /2 = cosh(£b). 
Correspondingly we have 
Q(0) 



l + M^/2 = cosh(£ £ ) 



(20) 



(21) 



(22) 



for the energies of momentum k eigenstates. 

Therefore, (M — Eq)/M — > holds in the quantum 
continuum limit M — ► independently of the spatial lat- 
tice size N. Restoration of the dispersion relation in the 
limit M -> requires N(M) -> oo for M -4-0. 

For comparison with U(l) in the next section we use 
Ek and draw in Fig. [I] 2^ as function of N for Eq = +0, 
which should be interpreted as a very small Eq > corre- 
sponding via (|2"Tj) to a similarly small M > 0. These are 
the first and the last of the curves shown. Respectively, 
the left and the right ordinate apply, the latter enlarging 
the approach to zero. 

With Eq = 1.0 and Eq = 2.0 two more examples 
are given in Fig. [TJ The left ordinate applies and the 
curves demonstrate the counter-intuitive result that the 
finite size corrections increase with decreasing correla- 
tion length £ = M~ x : £ = 1.04 for E = 1.0 and 
£ = 0.567 for Eq = 2.0. The approach to zero relies 
on k = 2tt/N -> 0. If a limit is considered for which |fc| 
stays finite for N — !> oo, the lattice dispersion is never 
fulfilled for energies calculated from the two-point corre- 
lation functions. This is obvious, because the lattice size 



enters into the definition (|T8|) of Meg only through an 
eventual dependence of the momenta on N. 



III. U(l) LATTICE GAUGE THEORY 

We consider U(l) LGT with the Wilson action Q on a 
4D hypercubic lattice with periodic boundary conditions 



S({U}) = J2Sa 



(23) 



with S a = ^e{U lin U ni2 U i2n Uj 2 i l ) 1 where ii, ji, i 2 and 
j2 label the sites circulating about the square □ and 
the Uij are complex numbers on the unit circle, t/y = 
exp{i(j>ij), < 4>ij < 2tt, Uji = . Expectation values 
are calculated with respect to the partition 



Z = 



! I <lo, exp[-PS({U})\ 



(24) 



where the product is over all links I = ij with ij near- 
est neighbor sites of the lattice. Wilson concluded that 
at strong couplings (small j3) the theory confines static 
test charges due to an area law for the product of gauge 
matrices around closed paths (Wilson loops), whereas it 
describes a massless photon at weak coupling (large /3). 
In between there is a transition between the confined and 
the Coulomb phase for which the ft value has meanwhile 
been accurately determined by MCMC simulations, see 
for instance Q. 

However, it turns out that correlations of Wilson loops 
in zero-momentum representations of the cubic group do 
not give any convincing signal for the existence of a mass- 
less photon. The presumed reason is a noisy power law 
fall-off, which cannot be followed beyond one or two steps 
in the lattice spacing a. A remedy found in Q relies on 
the use of the dispersion mass (U]) for the representa- 
tion of the cubic group. As reviewed in our introduction, 

ntlyj 
with zero [3|, £ 

In view of the deviation (|10l) from zero we decided 
to follow the finite lattice size dependence of (E2 is ) 2 
from small 4 3 16 to large 32 3 128 lattices. Using the 
Metropolis-heatbath algorithm of [l(| we performed 
MCMC calculations for which more details can be found 
in Q and measured correlation functions of the plaquette 
operator with momentum k\ ([5]) and in the repre- 
sentation of the cubic group. 

Our results at (3 = 1.05, 1.1 and 1.5 are summarized 
in Fig. [2j The left ordinate applies to the upper three 
curves. Using in Eq. (HJ) Eq = and {El) 2 as defined 
after Eq. ([5]) they show (Ej is ) 2 on a scale that includes 
the estimates from the smallest N = 4 lattice. With in- 
creasing lattice size an approach to zero is observed. The 
lower three curves, to which the right ordinate applies, 
reveal details. The lines are fits to a behavior of the data 
assumed to be governed by two exponential corrections 



until recenthylTj estimates of (Ej is ) 2 appeared consistent 
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FIG. 2: Lattice size dependence of the dispersion energy (|4|) 
for the photon in U(l) LGT. 



TABLE I: Estimates of the fit parameters (|25f) (error bars 
are given in parenthesis and apply to the last digits) and the 
goodness of fit Q. 



p 


ai 


a 2 


<Z3 


04 


Q 


1.05 


-1.84 (72) 


0.615 (83) 


0.0117 (19) 


0.0128 (60) 


0.49 


1. 10 


-1.91 (53) 


0.556 (54) 


0.0035 (19) 


0.0130 (20) 


0.31 


1.50 


-1.96 (46) 


0.547 (44) 


0.0002 (16) 


0.010 (37) 


0.89 



rized in Table Q] Large error bars are possible, because 
the fit parameter values are correlated with one another. 

The enlargements on the the right side reveal for the 
lower two /? values a considerable overshooting of the 
(-^dis) 2 = value, so that the final approach to zero is 
from above. Within our numerical precision it is unclear 
whether this overshooting continues for large f3 values. 
For f3 = 1.5 the coefficient 03 listed in the table is still 
positive, but well consistent with zero. 



IV. SUMMARY 

The negative value (fT0|) found in for the U(l) dis- 
persion energy squared is a lattice regularization effect, 
which we have derived with Eq. ([22]) analytically for a 
free scalar field. For Eq = and k = 2ir/N the disper- 
sion energy E^ is Q approaches zero with increasing lat- 
tice size for the free scalar field as well as for U(l) LGT. 
A distinction is the overshooting of the zero value, which 
is seen in Fig. [2] for U(l) LGT, but not in Fig. [Q for the 
free field. An astonishing feature of the free field exposed 
in Fig. [T] is the increase of the finite lattice size correc- 
tion with decreasing correlation length. The dispersion 
relation can only be restored for fcj — > 0, which includes 
the quantum continuum limit — ¥ and ki — > with 
~E % and k, t of Eq. © fixed. 



to zero 



( E Ii s ) = 01 exp(-a 2 N) + a 3 exp(-a 4 N) . (25) 
The fit parameters and the goodness of fit Q are summa- 
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